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■ We present a semiclassical theory of weak disorder effects in small struc- 

■ tures and apply it to the magnetic response of non-interacting electrons con 
Q\ ' 

fined in integrable geometries. We discuss the various averaging procedures 
^ describing different experimental situations in terms of one- and two-particle 



Green functions. We demonstrate that the anomalously large zero- field sus- 
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O ■ ceptibility characteristic of clean integrable structures is only weakly sup- 



pressed by disorder. This damping depends on the ratio of the typical size of 



the structure with the two characteristic length scales describing the disorder 
(elastic mean-free-path and correlation length of the potential) in a power- 
law form for the experimentally relevant parameter region. We establish the 
comparison with the available experimental data and we extend the study of 
the interplay between disorder and integrability to finite magnetic fields. 
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I. INTRODUCTION 



Electronic mesoscopic systems offer nowadays the possibility of being used as a labora- 
tory for studying quantum chaos. The main question of this novel discipline - the quantum 
signatures of the underlying classical dynamics - can be addressed in microstructures defined 
on high mobility semiconductor heterojunctions. This connection presents a considerable 
challenge to experimentalists since it implies complicated fabrication processes and delicate 
measurements. The challenge for theoreticians is not lesser since semiconductor microstruc- 
tures are very rich condensed matter systems (involving effects of temperature, confinement, 
disorder, electron-electron and electron-phonon interactions, etc) where the applicability and 
validity of simple models has to be clearly established. 

Within the simple model of a particle-in-a-billiard, important differences have been 
predicted^, and later measurediJl, in the transport through chaotic and integrable geome- 
tries. In the former nearby trajectories diverge exponentially and periodic orbits are usually 
isolated; the latter are characterized by having as many constants of motion in involution as 
degrees of freedom, and periodic orbits are organized in families on invariant torii. Chaotic 
cavities exhibit a universal behavior for the conductance fluctuations and weak-localization, 
characterized by a single scale. On the contrary, integrable cavities do not show generic be- 
havior presenting more fine-scale fluctuations and a non-lorentzian line-shape of the low-field 
magneto resistance. In the case of thermodynamical properties like the magnetic suscep- 
tibility, the differences between chaotic and integrable billiards are more spectacular since 
they involve an order-of-magnitude enhancement of the low-field susceptibility of integrable 
geometries compared to that of chaotic one Unlike the transport problem, the predicted 
different behavior according to the integrability of the underlying classical mechanics has 
not been experimentally confirmed yet. 

The residual disorder present in actual microstructures plays a special role in the quantum 
chaos studies. Indeed, any perturbing potential, such as the one provided by the disorder, 
immediately breaks the integrable character of the classical dynamics. Since small amounts 
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of disorder are unavoidable in actual microstructures, the question whether or not integrable 
behavior should be observed, naturally arises. It is then of foremost importance to establish 
if the differences between chaotic and integrable geometries persist when we go beyond 
the particle-in-a-box model. This interplay between integrability and disorder is the main 
subject of this paper. 

We start by characterizing the disorder. One limiting case is the absence of it, where 
the dynamics is determined by the non-random confinement potential (particle-in-a-box or 
clean models). On the other extreme we have the diffusive limit where the electron motion 
is a random walk between the impurities and the confining effects are not important. The 
strength of the disorder in the diffusive case is characterized by the transport mean free 
path l T : the mean distance over which the electron momentum is randomized. When l T 
becomes of the order of the typical size a of the microstructure, confinement and disorder 
are relevant. For l T > a we arrive at the ballistic regime where electrons can traverse the 
structure with a small drift in their momentum (going along almost straight lines), and 
their dynamics is mainly given by the bounces off the walls of the confining potential. In 
the ballistic regime the underlying classical mechanics still depends on the geometry and we 
would like to understand the different role of disorder in integrable and chaotic geometries. 

For short range impurity potentials (as typically found in metallic samples) the scat- 
tering is isotropic (s-type) and the momentum is randomized after each collision with an 
impurity. There is therefore only one length scale, namely It, characterizing the disorder. 
For smooth impurity potentials (as typically realized in high-mobility microstructures) the 
scattering is forward directed and l T may be significantly larger than the elastic mean free 
path / associated to the total amplitude diffracted by the disorderi. The regime l T > a > I 
is particularly interesting because it is ballistic (since the classical mechanics is hardly af- 
fected by disorder), but the single particle eigenstates are short lived. In a more technical 
language that we will precise on the sequel, we have I as given by a single-particle Green 
function and l T as by a two-particle Green functioni. We will study the interplay between 
disorder and confinement for physical observables that depend on one- and two-particle 
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Green functions, concentrating on the magnetic susceptibility of individual and ensembles 
of ballistic microstructures. 

The natural tools to attack the interplay between disorder and confinement are semi- 
classical expansions since they transparently convey at the quantum level the information 
about the classical mechanics. Supersymmetry0 and random matrix theories are quite pow- 
erful methods that have been widely used in recent studies of quantum chaos and disordered 
systems0H3, but are not applicable to our regime of interest since they deal with the ergodic 
universal (long time) properties of completely chaotic systems. Diagrammatic perturbation 
theory for the disorder can describe the diffusive regime^!, but calculations become exceed- 
ingly complicated when the confinement and the detailed nature of the impurity potential 
has to be considered. 

In our semiclassical approach we emphasize the dependence of disorder effects on the ratio 
between the finite system size a and the disorder correlation length £, showing that confined 
systems exhibit strong deviations from the bulk-behavior. In particular we demonstrate that 
for integrable geometries the effect of smooth disorder results in a power-law damping of the 
two-particle Green function properties, and we compare this behavior with that expected in 
chaotic systems. For completeness of the presentation we first briefly review in Sec. |TJ our 
work on the magnetic response of clean systemsiil. We then develop in detail a treatment 
of disorder in ballistic microstructures extending some preliminary work@. In Sec. |T| we 



present the disorder model and some general implications at the level of one- and two- 
particle Green functions. In Sec. [IV| and [V] we specialize in the impurity averaged magnetic 



susceptibility for individual and ensembles of microstructures. 



II. ORBITAL MAGNETISM IN CLEAN SYSTEMS: A BRIEF REVIEW 
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A. Thermodynamic Formalism 



In this section we present the basic thermodynamical formalism for obtaining the orbital 
magnetism within a semiclassical approach. We indicate the main ideas for its application to 
clean microstructuresi'i which will be further developed in the Sections f[V]|V] in order to allow 
the treatment of static disorder. The principle is to derive thermodynamical expressions for 
the free energy and the grand potential using a semiclassical approximation for the density 
of states. This allows us to calculate physical observables such as the magnetic susceptibility 
for the canonical and grand canonical ensembles. 

For a system of electrons confined to an area A at temperature T and subject to a 
perpendicular magnetic field H, the free energy F(T, H, N) for a fixed number N of electrons 
and the grand potential Q(T,H,n) (representing the coupling to a particle reservoir with 
chemical potential //) are related by means of the Legendre transform 

F(T, H, N) = fiN + Q(T, H, //) . (2.1) 

The canonical (x) and grand canonical (x GC ) susceptibilities of a confined electron gas are 
given by 

,= -i(4» • x =-44/ (2 - 2) 

The grand potential can be expressed in the form 

Q(T, H, fj) = -~ J dE d{E) ln[l + exp(/?(// - E))\ (2.3) 

(with (3 = l/k B T) in terms of the single-particle density of states d(E) which we decompose 
into a smooth mean and oscillating part according to 

d(E) = d{E) + d osc {E) . (2.4) 

As has been first noticed in the context of persistent currents in disordered ring S a dis- 
tinction between \ and x GC may be of crucial importance in mesoscopic thermodynamics: 
Although the number of electrons can be large for a mesoscopic system, the fact that N is 



fixed must be taken into account (by working in the canonical formalism) if a disorder or 
energy averaged magnetic response of an ensemble of isolated micro-systems is examined. 
According to Imry0 a convenient representation for the canonical free energy in terms of 
grand canonical quantities is obtained by expanding the relationship ( |2.1| ) to second order 
in \l — p, with a mean chemical potential p being implicitly defined by accommodating N 
charge carriers to the mean number of states 

N = N(p) = N(p) . (2.5) 

Here 

poo 

N(fi) = / dE d(E) f(E-fi) (2.6) 
Jo 

with the Fermi distribution function 

/(fl -^ i + e* Pl ,^-M)r (27) 



N is obtained in Eq. (2J3) by replacing d(E) by d(E). This finally allows an expansion of 
the free energy as^ 

F(N) ~ F° + AF (1) + AF (2) , (2.8) 

with 

F° = pN + n(fi) , (2.9a) 
AF (i) = fi osc^ ^ ( 2 Qh ^ 

AF(2) = Wp7) ( N ° sc ^ 2 ■ ( 2 - 9c ) 

The functions Q osc (p) and N osc (p) are expressed by means of Eqs. (|2.3| ) and ( p.6|) , respec- 
tively, upon inserting the oscillating part d osc (E) of the density of states (|2.4j) . The leading 
order contribution to F is given by the first two terms F° + AF^ yielding the susceptibility 
calculated in the grand canonical case. F° gives rise to the (two-dimensional) diamagnetic 
Landau-susceptibility which for billiard like systems is expressed as for the bulk as 



with g s = 2 the spin degeneracy. 



B. Semiclassical treatment of susceptibilities 

For a semiclassical computation of AF^ and AF^ and their derivatives with respect 
to H we calculate d osc (E, H) from the trace 

d(E, H) = --Im / (1tG e (t, r) (2.11) 

7T J 

of the semiclassical one-particle Green function. Its contribution to d osc (E) is given byi 

G s (r',r) = ^Aexp[^|-^0] , (2.12) 
as the sum over all classical paths t (of non-zero length) joining r to r' at energy E. 

S t = [ pdq (2.13) 



c, 



is the classical action integral along the path Ct- The amplitude D t takes care of the classical 
probability conservation, and r\t is the Maslov index. 

The evaluation of the trace integral ( |2.11| ) for chaotic and integrable systems leads to 



the Gutzwiller-@ and Berry-Tab or-El periodic-orbit trace formulas, respectively. In order to 
calculate the magnetic susceptibility at small fields one has to carefully distinguish! between 
the three possibilities of a chaotic billiard, the special case of an integrable billiard remaining 
integrable upon inclusion of the if -field, and the more general case where the field acts as 
a perturbation breaking the integrability of a regular structure. 

Since our main interest in the Sections p|JIV| and |V| will be devoted to disorder effects on 



the susceptibility of billiards being integrable at zero if -field we will focus here on the last 
case. There neither Gutzwiller- nor Berry-Tabor-trace formulas are directly applicable and, 
following Ozorio de Almeida^!, a uniform treatment of the perturbing if -field is necessary. 
In the integrable zero-field limit each closed trajectory belongs to a torus I-m and we can 
replace r in the trace integral Q2.11| ) by angle coordinates 0i specifying the trajectory within 
the (one-parameter) family and by the position 02 on the trajectory. For small magnetic 
field the classical orbits can be treated as essentially unaffected while the field acts merely 
on the phases in the Green function in terms of the magnetic flux through the area Am(®i) 



enclosed by each orbit of family M. Evaluating the trace integral ( |2.11|) along B2 for 
the semiclassical Green function of an integrable system leads in this approximation to a 
factorization of the density of states 

d° sc (E) = £ C M (H)dl(E) (2.14) 

into the contribution from the integrable zero-field limit 

d° M (E) = D M cos (k F L m - r/ M | - ^) (2.15) 

(Lm is the length of the orbits of family M and .Dm the semiclassical weig ht§) and the 
function 

1 f 27T 



fZTT 

l (H) = — d0i cos 
Ztt jo 



ZTT- 



(2.16) 



$0 

containing the if -field dependence ($0 = hc/e). Calculating AF^ from Eq. ( |2.9b| ) and 



taking the derivatives with respect to H gives the grand canonical contribution to the sus- 
ceptibility at small magnetic field 



Here, tm is the period of a closed orbit of family M and 



/ s t/t c h3 , . 

Mr) = . ' , r c = J± 2.18 

smh(r/r c ) 7r 



is a temperature damping factor which arises from the convolution integral in Eq. ( |2.3| ) and 
gives an exponential suppression of long orbits. This is important from a physical as well 
as computational point of view, as conceptual difficulties associated with the questions of 
absolute convergence of semiclassical expansions at zero temperature do not arise. 

Eq. ( p,17| ) is the basic equation for the susceptibility of an individual microstructure. 



When considering ensembles of ballistic microstructures however, an average ( — " ) over 
energy (i.e. k F ) or over the system size a has usually to be performed and leads to variations 
in the phases (actions S/h = /c f Lm) of the density of states (|2.15|) which are much larger 
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than 27r. Therefore, x vanishes upon ensemble average. In order to characterize the orbital 
magnetism of ensembles we introduce the typical susceptibility x (t) = (x 2 ) 1//2 (the width of 
the distribution) and the ensemble average x (its mean value, which is non-zero because 
of the term AF^ in the expansion Eq. Q2.8|) ). They are of theoretical interest because of 
being based on two-particle Green functions and they are relevant for the description of 
experiments on ensembles of mesoscopic systems. 

If we assume that there are no degeneracies in the lengths of orbits from different families 
M we obtain for x {t) 

In calculating x, the grand canonical contribution x^ from AF^ vanishes under energy 
average and the canonical correction AF^ in Eq. fl2.8|) gives in semiclassical approximation 
using Eq. ( p.9c|) 



Xl Xl gs 2 \2nA 



E 1 ^^) 2 ^- (2-20) 

The Eqs. fl2.17|) -( p.20|) provide the general starting point for a computation of the suscepti- 
bility of integrable billiards at small fields. 

As an important example, which is also of experimental relevance0, we will apply the 
results to square billiards. At finite temperature x is essentially given by the family M = 
(1,1) of the shortest, flux-enclosing periodic orbits depicted in Fig. |I[ A complete treatment 
including families of longer orbits is given in Ref.i. Instead of Gi we use the lower reflection 
point xo as orbit parameterization within the family. The orbits (1,1) have the unique length 
Lxi = 2y / 2a and enclose a normalized area A(xq) = 4nxo(a — Xo)/a 2 . Computation of d^i(fJ>) 
for the square geometry gives for x^ (Eq. ( |2.171 )) 

f a dxc 



X° 



£ ^A\x Q ) cos(M(*o)) sm (k r Lu + \ ) (2-21) 



as a function of the total flux ip = Ha 2 /§ with $ — hc/e. The prefactor 

x ° = x '(7l)^ (M)3/2flT(ill) ' (2 ' 22) 
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shows the (fc F a) 3/ ' 2 -dependeiice typical for (nearly-) integrable systems. 

For the square geometry the Equations fl2.19|) and (|2.20 ) for the susceptibilities x (t) and \ 
(characterizing different ensemble averages) can be reduced to (including only the dominant 
contributions from the family (1,1)) 



X_ 

X G 



(t) 



X 



(1)' 



1 r d^o^ 2 cos((pA(x )) 
V2 Jo a 



and 



X_ 

x° 



1 r a dx 

2 Jo 



a dx' n 



(2.23) 



with 



a Jo a 



A 2 _ cos(yx4_) + A 2 + cos(</9^4_| 



^2.24) 



Xl (V2- 



7T 



(k F a) -R|(L n ) 



(2.25) 



and ^4± = A(xo)±A(x' ). Although the integrals ( |2.21| ),( pl23"D and ( p.24| ) can be evaluated 
analytically in the clean case (leading to Fresnel functions of the magnetic flux!) the above 
expressions serve as suitable starting points for the discussion of disorder effects on ensembles 
of microstructures discussed in Sections |V| and M. 



III. SEMICLASSICAL APPROACH TO WEAK DISORDER 

Disorder is usually studied in terms of the ensemble average over impurity realizations, 
since it is a perturbation of a electrostatic potential whose detailed nature is unknown. 
Typically, quantum perturbation theory is followed by the average over the strengths and 
positions of the impurities. This approach is suited for macroscopic metallic samples (which 
are self-averaging) or ensembles of mesoscopic samples (where different samples present dif- 
ferent impurity configurations). The possibility of measuring a single disordered mesoscopic 
sample posses a conceptual difficulty since there is not an average process involved. When 
discussing the effect of disorder on the orbital magnetism of microstructures, it is therefore 
necessary to distinguish between the behavior of an individual sample and an ensembleil. 
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Moreover, we have to consider the cases where the Fermi energy and size of the mi- 
crostructures are kept fixed under impurity average and the cases where these parameters 
change with the different impurity realizations. These various averages, that will be thor- 
oughly discussed in the remainder of the paper, can be expressed in terms of the impurity 
average of one- and two-particle Green functions. Therefore we perform in this section a 
general treatment of disorder effects on the basis of semiclassical expansions of Green func- 
tions. The Green function formalism, which is useful for a wide range of physical problems, 
can be applied to thermodynamical quantities like the magnetic susceptibility (Sec. fV] and 
|V]) as well as to quantum transport problems. 



1. Disorder models 

Our basic assumptions for the treatment of disorder are the following: We study a 
spatially random potential V(r) characterized by a correlation function 

C(|r - r'|) = (V(r)V(r')) (3.1) 

with a typical correlation length £ and a mean disorder strength C° = C(0). We will make 
use of a Gaussian correlation 

C(|r-r'|) = C°exp(-^^) (3.2) 

which allows us to derive analytical expressions for disorder averages considered belowll. 
The disorder correlation function (|3.2| ) can be viewed as being generated by means of a 
realization i of a two-dimensional Gaussian disorder potential given by the sum 

of the potentials of N{ independent impurities located at points Rj with uniform probability 
on an area V. The strengths Uj obey (ujUj,) = u 2 5jj>. The disorder strength (as defined in 
Eq. Q) is 
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with rti = iVj/V. For £ —>■ this model yields the white noise case of 5-function scatterers 
V(r) = Ef'^r-R,-)- We will use the model of Gaussian disorder for some analytical 
calculations and for numerical quantum simulations. However, the general results expressed 
in terms of the correlation function C(\r — r'\) will be valid for any kind of disorder. 

As we will show, disorder effects depend on several length scales: the Fermi-wavelength 
A F of the electrons, the disorder correlation length £ and the size a of the microstructure. In 
the bulk case of an unconstrained two-dimensional electron gas (2DEG) we will distinguish 
between short range (£ < A F ) and finite range (£ > A F ) disorder potentials. In the case of a 
microstructure a third, long range regime for £ > a > A F has to be considered. The cleanest 
samples used in today experiments are in the finite range regime a > £ > aJhI. 



2. Single-particle Green function 

If we assume a microstructure with size a ^> A F (a condition which is always met in 
lithographically defined samples) and work in the finite range or long range regime, where 
the disorder potential is smooth on the scale of A F , a semiclassical treatment is well justi- 



fied. A natural starting point is the semiclassical expression Q2.12 ) for the single-particle 
Green function Ge{t', r) as a sum over the contributions from classical paths. The classical 
mechanics of trajectories with length L t <C l T (the transport mean free path) is essentially 
unaffected by disorder. Therefore the dominant effect on the Green function in Eq. ( |2.12| ) 
results from shifts in the semiclassical phases due to the modification of the actions while 
the amplitudes D t and topological indices rjt are nearly unchanged. The first-order approxi- 
mation to the classical action ( |2.13|) along a path C t in a system with weak disorder potential 



is 



S d t ~ SI + 5S t , (3.5) 
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where the clean action Sf is obtained by integrating along the unperturbed trajectory C 4 C 
without disorder (i.e. S% = k F L t in the case of billiards without magnetic field) instead of 



the actual path C t . The correction term SSt is obtained, after expanding p = y2m[E — V(q)] 
for small V/E, by the integral 



SS t 



V(q) dq . 



(3.6) 



In this approximation an impurity average (. . .) acts only on 5S t and the disorder averaged 
Green function reads 



(G E (r',v)) = Y,G% t (r',r) (exp 



> ■ 



(3.7) 



Here G C E t is the contribution of the trajectory t to the zero-disorder Green function G C E . 

For trajectories of length L t ^> £ the contributions to 5S according to Eq. (|3.6| ) from 
the disorder potential at trajectory segments separated by a distance larger than £ are 
uncorrelated. The related stochastic accumulation of action along the path can be therefore 
interpreted as determined by a random-walk process, resulting in a Gaussian distribution 
of SS t (L t ). For larger £ or shorter trajectories (L t ^> £), one can still think of a Gaussian 
distribution of the de-phasing SS t provided V(r) is generated by a sum of a large number 
of independent impurity potentials. As a consequence of the Gaussian character of the 
distribution of 5S t (L t ), the characteristic function involved in Eq. ( |3.7|) is given by 



(exp 



n 



) = exp 



2tr 



(3.8) 



and therefore entirely specified by the variance 

(^ 2 ) = -i / dq [ dq'(V(q)V(cO) , (3.9) 

v F Jc c t Jc$ 

which is expressed as the mean of the disorder correlation function C(|q — q'|) when the 
unperturbed orbit is traversed. 

If we consider, to start with, an unconstrained 2DEG the sum in Eq. ( |3.7| ) is reduced to 
the direct trajectory joining r and r'. If L = |r — r'| ^ £ the inner integral in Eq. 
be extended to infinity and we obtain 



can 
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(5S 2 ) = ^ fdqC(q). (3.10) 

v f J 

The semiclassical average Green function for the bulk exhibits therefore an exponential 
behaviorlll'il (on a length scale l T > L ^> £) 

(G E (v',v)) = G c E (r',r) exp(-| ), (3.11) 

with the damping governed by an inverse elastic mean free path 

'dgC(q). (3.12) 



I k 2 v 2 P 

In the case of Gaussian correlation C(q) is given by Eq. fl3.2|) and we get 



l = JL^. (3.13) 



Using the disorder strength ( |3.4|) we have 



; = i^*. (3.M) 



In Appendix [A] we discuss the relation between the semiclassical elastic MFP's 
(Eqs. (|3.12|) - (|3.14|) ) and the MFP obtained from quantum diagrammatic perturbation the- 
ory for the bulk for the disorder model ( |3.3| ). The semiclassical and the quantum result 
(Eq. (POJ)) agree asymptotically to leading order in kp£. In the limit of small £, especially 
£ < Xp, our semiclassical approach is no longer applicable^. However, Eq. ( |3.11| ) still holds, 
but with I replaced by Is given in Eq. ( ]A4| ) . 

We now turn from the semiclassical treatment of the bulk to that of a confined system. 
In the constrained case in the limit l T a impurity scattering is the dominant process^. 
This gives rise to diffusive motion, and thus there is no essential difference to the bulk for 
the damping of the Green function. We will treat the ballistic regime l T >a where both, the 
confinement and the impurities have to be considered. A treatment of l T in the appendix 
shows that for finite £ the transport MFP I is considerably larger than the elastic one and 
a ballistic treatment is therefore well justified, even if I is of the order of the system size. 
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In contrast to the bulk case a disorder averaged confined system is no longer trans- 
lationally invariant and one has to impose in quantum calculations the correct boundary 
conditions of the geometry. Confinement implies semi-classically that G° E (r',r) is given as 
a sum over all direct and multiply reflected paths connecting r and r'; disorder modifies the 
corresponding actions according to Eq. (^.6|). 

In the regimes of short- and finite-range scatterers, the damping of each contribution 
(GrE,t) to (Ge) is given, analogous to the bulk expression (|3.11|) , (using Eq. (|3.10|) ) by 

(Ge(v',v)) = £G^(r',r) exp (-^) . (3.15) 

Here, L is now replaced by the trajectory length L t > a ^> £. This gives an individual 
damping exp(— L t /2l) for each geometry-affected path contributing to (Ge)- 

In the long range regime and for £ ~ a the correlation integral ( |3.9| ) can no longer be 
approximated (as for £ <C L t ) by L J*™ dqC(q) due to correlations across different sectors of 
an orbit (with distance smaller £). Therefore the orbit-geometry enters into the correlation 
integral. For (>awe can however expand C(\r — r'|) and obtain in the case of Gaussian 
disorder (up to first order in £~ 2 ) C(|r — r'| ~ C°[l — (r — r') 2 /(4£ 2 )]. In this approximation 
the integral ( |3.9| ) gives for the Green function damping an exponent 

(5Sf) 1 L, 2 / 1 I t \ , . 

■ (316) 

I t = {1/L t ) f c r 2 (q)dq can be regarded as the "moment of inertia" of the unperturbed 
trajectory C t with respect to its "center of mass" (1/L t ) J Ct r(q)dq. Eq. (gT6|) shows that 
the damping in the long range regime depends quadratically on L t (in contrast to linear 
behavior in the finite range case or bulk). The length scale of damping is now given by the 
geometrical mean of the bulk MFP I and £. The leading damping term does not depend 
on the specific orbit geometry since it essentially reflects the fluctuation in the mean of the 
(smooth) potentials of different impurity configurations. Inclusion of higher powers of 
leads to additional contributions from higher moments J c r n (q)dq on the RHS of Eq. (ft.l6Q . 
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3. Two-particle Green function 



Density correlation functions in general or the typical susceptibility (Eq. (|2.19|) ) and 
ensemble averaged susceptibility (Eq. (|2.20| )), which will be treated in the subsequent sec- 
tions, involve squares of the density of states. Writing the latter, Eq. ( |2.11| ), in terms of 
the difference between advanced and retarded Green functions (G + — G~) we are left with 
products of one-particle Green functions. The terms of most interest are the cross products 
G + (r,r') G~(r,r') = G + (r,r')G + *(r',r), because they survive the energy average and are 
sensitive to changes in the magnetic field. 

Since in the non-interacting approach we are using, the two-particle Green function 
factorizes into a product of one-particle Green function!! we will use the former as a synonym 
for the latter. The semiclassical average for products of single-particle Green functions will 
be quantitatively performed for the susceptibility of confined integrable systems in Section 
0, and we discuss here the underlying ideas for the general case. 

Considering for instance the product G(rx, r 2 )G*(r[, r' 2 ), the effect of the disorder poten- 
tial can be taken into account perturbatively for each realization of the disorder in the same 
way as before by Eqs. ( |3.5| )-( |31j| ). Using the same kind of argument, one can therefore write 
the disorder average as a double sum over the averaged contributions from trajectories t and 
t' 

(Ge G^>=EE <G*,t G* Et ,) = EE G% tr {e {i/mSt - 5St,) ) (3-17) 
t t' t v 

{(5S t -5S t ,) 2 ) 



EE G %,t g*,*«p 



2h 2 



It is necessary here however to take into account the correlation of the disorder potential 
between points on trajectories t and t'. One limiting case for instance would be that t and 
t' are either the same trajectory or the time reversal one of each other. In these cases their 
contribution acquires exactly the same phase shift and (Gej G* E t ) = \G c E j\ 2 . Within our 
approximation the diagonal contributions t = t', which e.g. are responsible for the classical 
part of the conductivity, remain thus disorder-unaffected, since we assume the trajectories 
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have a length much smaller than It- (A semiclassical consideration of these effects for 
trajectories of length of the order of It or larger was performed in Ref.H for the bulk, giving 
a damping of the two-point Green function on the scale of l T .) At the opposite extreme, 
if trajectories t, t' are completely uncorrelated, i.e., for long trajectories in classical chaotic 
systems or trajectories in integrable systems with a spatial distance larger than £, the average 
in Eq. (|3.17|) factorizes: (G^ G* Et ,) = (Gej) 4 (G Et >) and lead to single-particle damping 



behavior. 

The double sum Eq. (|3.17|) may however involve pairs of trajectories which stay within a 
distance of the order of £ (as for nearby paths on a torus of an integrable system) . In this case 
the behavior of (G E ,t G* Et ,) is more complicated and depends of the confinement geometry 
of the system under consideration. As a simple illustration of the interplay between disorder 
correlation and families of orbits, let us consider for the case of the bulk the product of 
G(ri,r2) joining r\ = (0,0) to r 2 = (L, 0) with G*^,^) joining r\ = (0,y) to r 2 = (L,y), 
with L ^> £ but y possibly of the order of £. Introducing the function 

r+oo 

K(y)= C(x,y)dx (3.18) 



(for Gaussian correlations Eq. (|3.2|), K(y)/K(0) = exp(— y 2 /(4£ 2 ))), the variance of the 
de-phasing is obtained as 

( W -^fl =8 L fflz£M (3 . 19) 

V 
<- - F 



and therefore (G E G* E ) = G% G c £f{y) with 



hv) = ex p 



L ( _ K(jj)\ 
l e { K(Q)J 



(3.20) 



The function f(y) expresses in a very simple way that as y — > 0, the effect of disorder 
disappears (/(0) = 1) while for y ^> £ the function f(y) behaves as the square of single 
particle Green function damping. 



17 



IV. FIXED SIZE IMPURITY AVERAGE OF THE MAGNETIC 

SUSCEPTIBILITY 



We consider here a disorder average (which will henceforth be called a fixed-size impu- 
rity average) of an ensemble of structures for which the parameters of the corresponding 
clean system (geometry, size, chemical potential) remain fixed under the change of impurity 
realizations. In Section [V], we will then treat the more realistic case of the orbital magnetic 
response of a combined energy (or size) and disorder average. 

As shown in the previous section, averages over weak disorder exponentially damp, but 
do not completely suppress oscillatory contributions (with phase k F L t ) to the single-particle 
Green function from geometrical paths in confined systems. An observable quantity depen- 
dent on these contributions is the disorder averaged susceptibility of an ensemble of billiards 
of the same size or same clean-system Fermi energy, which will be studied first. 

We will treat regular billiards at zero or small magnetic fields, where the integrability is 
approximately maintained and the density of states has the independence of the formulae 
( p,14j )- (|2.16| ). The general result for x > Eq. ( |2.17|) , formally persists with the replacement 
of C M by 

(C M (H)) = ^- r&Q x cos 

Z7T JO 

where (55^(Gi)) is given by Eq. ( p.9| ) with the integrals performed along the orbits of the 
family M parameterized by 0i. In the finite range case (if all orbits of a family M are of 
the same length as in billiards) each family exhibits a unique disorder damping giving a 
contribution 

(Xff) = xff-exp(-^|l) (4.2) 

to the ballistic susceptibility. Xm is the contribution of family M to the clean susceptibility 
(Eq. ( ETT7D ) and (SS M )/2h 2 = L M /2l. 

In the case of square billiards, where the dominant contribution stems from the family 
(1,1), we obtain, in analogy with Eq. ( |2.21 ), 
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exp 



2h 2 



(4.1) 



M ~ i^l = f ^±A 2 {x ) cos(^4(x )) /sin f fc F L 



7T | JS^o)^ 

11 



(4.3) 



4 h ) I 

with x° given by Eq. ( [2.22|) . For a square billiard f^S^o) is independent of x for the finite- 
as well as for the long-range regime since In = a 2 / 12 (entering into Eq. Q3.16Q ) is the same 
for all orbits (11). Therefore Eq. ( |4.2| ) with M = (1, 1) holds for both limiting cases. In the 
same way as for the damping of the one-particle Green function (Eq. ( |3.15| )) we obtain for 
square billiards at finite temperature in the finite range regime 

(x)^(x (1) ) = X?-exp(-^) , (4.4) 

where Xd denotes the susceptibility of the system without disorder. 

In order to control the validity of our analytical semiclassical approximations we 
performed numerical quantum calculations by diagonalizing the Hamiltonian for non- 
interacting particles in a square billiard subject to a uniform perpendicular magnetic field 
and a random disorder potential of the form of Eq. (|3.3| ). For a given selected correlation 
length £ and quantum mechanically calculated elastic MFP Z qm and fixed Fermi momen- 
tum k F the product of the mean number of impurities per area and squared mean impurity 
potential, n^w 2 , is determined by Eqs. (|A3| , |A4D . We found that our numerical results are 
essentially independent of the choice of rii (with u 2 adjusted accordingly) for rii > 200 and 
used this value for the calculations presented here. The positions Rj of the impurities were 
chosen as independently distributed and for the Uj we used a box distribution. 

Each impurity configuration a has a self-averaging effect for an individual square billiard 
(for £ < a) due to the differences of the impurity potential V a (r) across the structure. In 
an average over an ensemble of square billiards, differences in the mean impurity potential 
V a = (1/a 2 ) / dr V a (r) (the integral is taken over the area of the billiard) between different 
squares lead to an additional damping. It is characterized by the variance 



(V) 



2, U 2 Ui 



a?rf 
u 2 rij 



1 / _ 2 

f] erf (77) H — = (e v — 1 

7T V 



1 2 



4vr£ 2 
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a 

V=2£ (4-5) 
for £/a — > 00 (4.6) 




for f 







(4.7) 



a? 



In the limit of £ ^> a our numerical calculations showed that the self-averaging effect is 
negligible (since the impurity potential is essentially flat across the square) and the clean 
susceptibility of an individual structure remains practically unaffected by disorder. In this 
limit variations in the mean potential V of an ensemble (Eq. (|4.6| )) dominate the damping. 
In the limit of white noise disorder, fluctuations in the mean V of different samples play a 
minor role and self-averaging is the predominant process for an integrable system: In semi- 
classical terms different trajectories of a family of closed orbits are perturbed by white noise 
disorder in an uncorrelated manner. Therefore we do not observe considerable differences 
between the susceptibility of a single disordered billiard of integrable geometry and the 
corresponding ensemble for ( < a. In a chaotic billiard this self-averaging effect does not 
exist (for not too small £, see end of Sec. [V])), since orbits are isolated. Therefore distinct 
differences between an individual disordered sample and an ensemble of disordered billiards 
are expected. 

To improve the statistics of our numerical ensemble average for square billiards we per- 
formed an average over disorder configurations with the same mean V and in addition 
averaged over V according to Eq. (|4.5| )il. Fig. |2| shows results of the numerical quantum 
simulations for the average susceptibility (%) of an ensemble of squares with fixed size but 
different disorder realizations at a temperature k B T = 3g s A, where A is the mean level 
spacing. The characteristic oscillations in k F a show an interchange between para- and dia- 
magnetic behavior on a scale k F Lu. This indicates that they are dominated by contributions 
from the shortest flux-enclosing orbits of the family (1,1) (according to Eqs. ( |2.171 , [4.3| )), as 
has been already shown for the clean case in RefsH'i. Fig. § demonstrates the damping of 
the clean susceptibility (dotted line) with decreasing elastic MFP I /a = 4, 2, 1, 0.5 for fixed 
£/a = 0.1 (which represents a typical disorder correlation length in experimental realiza- 
tions). Variations in the mean V lead to a de-phasing of the oscillations in the finite range 



case on a scale (Sk)a ~ (4vr) 1/,4 w^/Z qm (^) which is, as discussed above, small compared to 
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the self-averaging effect in this regime. 

Fig. |3] depicts the quantitative comparison between numerical and analytical results: 
It shows the logarithm of (x) normalized to the corresponding zero-disorder susceptibility 
as a function of the inverse MFP for different correlation lengths £. The semi-classically 
predicted exponential damping (Eq. ( |4.2j )) is shown as straight lines for the short range 
(£ <C a, Eq. (|3.11|) , full line for £ = 0) and long range (£ > a, Eq. (|3.16|) , dotted lines 
for £/a =4,2,1 from the top). The semiclassical predictions accurately agree with the cor- 
responding quantum results (symbols) for £/a = 4,2,1,0 and fail for intermediate values 
£/a = 0.5,0.2 (squares and diamonds) which are off the range of validity of the approxi- 
mations. The transition from self-averaging dominated (£ — > 0) suppression to damping 
according to fluctuations in the floor V (for £/a — > oo) turns out to be non-monotonic. 



V. COMBINED IMPURITY AND ENERGY AVERAGE OF THE 

SUSCEPTIBILITY 

In currently experimentally realizable structures disorder averages cannot be performed 
independently from size-averages since the detailed features of the confining potential do not 
remain unchanged for different impurity configurations. From the basic expressions Q2.21|) 



and (|4.3| ) for the susceptibility we see that changes in size a give rise to rapid variations in the 
phase k F a (on a quantum scale) and a much slower secular variation through the geometrical 
factors A. Thus, the effect of small size variations is equivalent to an energy (k F ) average. 
As discussed in Sec. |IJ for the clean case, variations in k F lead to vanishing x ■ Therefore 
we have to use the typical and energy averaged susceptibilities (see Eqs. ( |2.19| ) and (|2.20|) 



for their definition in the clean case). When disorder is introduced we must consider energy- 
and disorder averages. The typical susceptibility is now defined by = (x 2 ) 1 ^ 2 - It applies 
to the case of repeated measurements on a given microstructure when different impurity 
realizations (and simultaneous changes in k F ) are obtained by some kind of perturbation 
(e.g. cycling to room temperature). From now on we will reserve the term xti f° r the clean 
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typical susceptibility (x 2 ) 1 ■ The energy and impurity averaged susceptibility (%) describes 
the magnetic response of an ensemble of a large number of microstructures with different 
impurity realizations and variations in size. This is the situation of the experiment of Ref.1^11 
that we discuss in the sequel. 



A. Integrable systems: The square billiard 

The semiclassical results for x® and (x) for a system of integrable geometry are obtained 
in an analogous way as we proceeded for (x) in Section [IV], that is by including in the integral 
fl2.16| ) for Cm a 6i-dependent disorder-induced phase exp(i5S(Qi)/h) (see Eq. (|4~l|)). How- 



ever, now we have to take the square of C M (respectively d 2 C-M./dH 2 ) before the impurity 
average and cross correlations between different paths B and 0' on a torus M or between 
different tori have to be considered. We discuss this effect, typical of integrable systems, for 
the case of a square billiard. For sake of clarity we moreover assume a temperature range 
such that only the contribution of the shortest closed orbit has to be taken into account. 
Instead of Eq. (|2.23|) and ( p.24|) which hold for the clean case, the contribution of orbits of 
topology M = (1, 1) for the typical susceptibility now reads: 
/ \ 2 1 f a dx 1 ' 



\X° J 2 Jo a Jo a 
with x° defined as in Eq. ( |2.22| ). The function 



£ A 2 {x )A 2 {x' ) cos (<pA(x )) cos (<pA(x' )) f(x , x' ) , (5.1) 



/(*„, A) = ^exp (6S(x ) - (5-2) 



exp 



1 (5S 2 (x )) + (5S 2 (x' )) - 2(5S(x )5S(x' )) \ (5.3) 



2h 2 

accounts for the effect of disorder on pairs of orbits x and x' . See Eq. ( |3.20| ) for the 
treatment in the general case. For the magnetic response of an energy- and disorder- 
averaged ensemble we find correspondingly: 

/ — / ^ U 2 _cos(<^)+*4 2 + cos(<M + ) ]f(xo,x'o) (5-4) 
Jo a Jo a 1 J 

with x° defined in Eq. ( ggg ) and A± as in Eq. (^24|) . 
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(x) 1 f a dx f a dx' 
Y = 2 



1. Short range case 



We begin with the discussion of the short range case: Although we reach the border 
of applicability of our semiclassical approximation for £ — > 0, it shows us that in this 
limit orbits with xq ^ x' are disorder-uncorrelated and all such pair contributions are 
exponentially damped. Using exclusively the family (1,1), one obtains an overall suppression 
of the typical and average susceptibility at finite temperature according to 

limx (t) =Xci e- Lll/2ls (5.5) 
lim(x) = x e~ Lll/ls . (5.6) 



Note that the exponent for (%) differs by a factor 1/2 from that for (%) (see Eq. ( |4.2| ) and 
subsequent text). 

Fig. [| depicts the k F a dependence of the ensemble averaged susceptibility (x) in the 
short range case £ = 0. The dotted curves showing the semiclassical analytical formula ( |5.6| ) 
are compared with a direct quantum mechanical calculation of (x ) (using the numerically 
obtained N oso (jl) in Eq. (|2.9cf )) for disorder ensembles of different impurity strength equiv- 



alent to an elastic MFP l$/a = 00, 8, 4, 2 and 1 at k F a ~ 65 (from the top). Note, that the 
effective MFP decreases along the curves with decreasing k F (see Eq. ( |A4|) ) and the localized 
regime may eventually be reached for small k F a. At the limit of the ballistic regime at small 
I ~ a the semiclassical result begins to differ from the quantum one although the functional 
behavior remains the same. This arising difference may be related to non-ballistic scattering 
from impurities which is not included here. 



2. Finite range case 

In the finite range A F < ( <C a, the phase shifts SS(x ) and SS(x' ) in f(x ,x' ) are 
accumulated in a correlated way, if the spatial distance of two orbits xq and x' Q is smaller 
than £. To evaluate the product term 2{SS(x )SS(xq)) in the exponent of f(x ,x' ) in this 
regime the integrations are performed as in Eq. ( |3.9|) but with q and q' running along 
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paths starting at xq, respectively x' . Ignoring the additional correlations occuring near the 
bounces off the boundaries of the billiard, the trajectories x and x' Q (see Fig. |l|) can be 
regarded as straight lines remaining at a constant distance y — \x — x' \/^/2 from another. 
We can therefore approximate f(xo,x' ) by f(\xo — x' \/\^2) with the function / given by 
Eq. ( |3.20D . For Gaussian correlation we thus have 



f(x ,x' )=ex P > ~" ' ' (Xo X>o)2 



Ln 
I 



cxp 



(5.7) 



Orbits separated by \x — x' \ £ are disorder-uncorrelated and exponentially suppressed: 
J(xq, x'q) ~ exp(— L\i/l). For those orbits the individual random disorder leads to an uncor- 
related detuning of the phases. In contrast to that, disorder only weakly affects trajectories 
separated by \xq — x' \ < £. 

The disorder averages in the finite range regime lead, by means of the function /, to a 
non-exponential damping of the susceptibilities for systems with families of periodic orbits. 
This behavior becomes obvious for the case of square billiards where at H = the integrals 
( |5.1D and ( |5.4|) can be evaluated analytically in the limits of Ln I (extreme ballistic) and 



L\i 3> / (deep ballistic). We find for the typical and average susceptibility at H = in the 
finite range case for Ln <^ I 



(*) / — / 



1 - 1 - ct - , (5.8a) 



a 



X I \ a 

and for Ln 3> I (by steepest descent): 



<X) - Ln( l-cJ-), (5.8b) 



f-(!)(0 /2 - 



The constants in the above equations are c t = (20/7) v27r and c a = 2v2vr. Eqs. ( p.8| ) express 
the limit of very weak disorder, showing that the small disorder effect is further reduced 
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due to the correlation of the disorder potential. The other limit, Eqs. (|5.9|) , is noticeably 
more interesting since it shows that disorder correlation effects lead to a replacement of the 
exponential disorder damping by a power law. 

Fig. [| depicts in logarithmic representation our collected results for the disorder averaged 
typical (a) and averaged (b) susceptibility for square billiards (at H = and k B T = 2g s A) as 
a function of the inverse elastic MFP for different disorder correlation lengths. The symbols 
denote results from numerical quantum simulations described in the previous section and the 
full curves semiclassical results from numerical integration of the Eqs. (|5.1|) and (|5.4j ). For 
the short range case £ = they reduce to Eq. (|5.6|) predicting an exponential decrease with 
exponent Ln/l which is in line with the quantum calculations (circles). The semiclassical 
results for the finite range are on the whole in agreement with the numerical results for 
£/a = 0.1 (diamonds), £/a = 0.2 (triangles) and £/a = 0.5 (squares). The semiclassical 
curves seem to overestimate the damping of the typical susceptibility. The dotted curves 
(shown for a/l > 1) depict the analytical expressions ( |5.9| ) in the regime L n > I. Since 
for finite £ the transport MFP It > I (see Eq. ( |A6| ) ) , this regime can still be considered as 
(deep) ballistic and our semiclassical assumptions being based on straight-line trajectories 
remain valid. 

As the semiclassical formulae already indicate, the overall disorder behavior of (x) and 
X® is quite similar. 

3. Long range case 

For completeness, we will consider the effect of the disorder for the long range regime: 
We can use the Eqs. ( |5 . 1| ) and (|5.4j ) but cannot calculate the disorder function f(xo,x'o) 
in the same way as for the finite range. We can however, similar as for (x) in Section [TV] , 
expand the exponent —((5S(xq) — 5S(x' )) 2 ) of f(xo,x' ) in Eq. ( |5.2| ) for small a/£. In the 
case of the square all orders up to (a/£) 8 vanish and we find a very small overall reduction 
of the clean averaged susceptibilities (from family (11)) given by 
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For square billiards this leading order contribution does not depend any longer on Xq. The 
energy- and disorder-average (x) exhibits the same damping as (x^) 2 - Note that besides 
the high order in (a/£) the prefactor is rather small. This weak suppression of the averaged 
susceptibilities can be related to the fact that in the long range case, different sectors of 
the contributing periodic orbits are highly correlated. As visible in Fig. ||(a), the quantum 
mechanical results (squares) for x at £/a = 0.5, which are closest to the long range case, 
exhibit already a very weak damping. 

B. Disorder effects at finite // field: from integrable to chaotic behavior 

In Fig. |] we compare the ratio {x /x°Y (obtained from calculating the integral in 
Eq. ( |5.1D ) as a function of the dimensionless flux ip = Ha 2 /$> for the clean case and for 
disorder characterized by I = a and £ = 0.1. This figure shows that the damping due to 
disorder is maximal at zero field, but that already for (p = 5 the disorder seems not to affect 
the magnetic response any further. 

The origin of this behavior can be understood readily by observing that as soon as <p 
is larger than one, the integral Eq. (|5.1|) is correctly approximated by a stationary phase 
approximation!. The stationary point x s Q = a/2 corresponds to the two periodic orbits of 
the perturbed system, and only the trajectories such that 

(x - x s ) 2 ip < 1 (5-11) 

actually contribute to the integral. The magnetic field causes a de-phasing of the con- 
tributions of the various trajectories of the family, thus breaking the integrability of the 
system. This effect is responsible for the overall decrease of the typical susceptibility as 
the field increases. In this respect clean and disordered square billiards are not equivalent. 
In the disordered case, trajectories separated by a distance larger than £ are already not 
contributing in phase. Therefore the additional magnetic field affects the magnitude of the 
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susceptibility much less. This remains true up to the point where the condition ( |5.11 



im- 



plies |xo — Xq\ < £ in which case the disorder is not effective anymore, and the two curves 
coincide. 

Therefore the behavior of the disorder damping we discussed in the previous subsection 
is characteristic for integrable geometries. For chaotic systems diagonal contributions (pair 
products of the same periodic orbit) are barely affected by disorder. This behavior is similar 
to that of an integrable systems at finite field. When evaluating the contribution to the trace 
of the Green function in the neighborhood of a periodic orbit by stationary phase approx- 
imation, (as for the derivation of the Gutzwiller trace formula) only orbits extremely close 
to the periodic orbit under consideration actually contribute. Unless £ is taken exceedingly 
small, all these trajectories will see the same disorder potential. 

As a final remark, note that non-diagonal contributions (pairs of different paths) are fully 
damped upon impurity average for chaotic as well as integrable systems, since the disorder 
potential along two different trajectories is usually not correlated (see also the discussion 
of the averaged Green function product after Eq. ( |3.17| )). Therefore non-diagonal contri- 
butions, which might be necessary to be considered in the clean case@, are exponentially 
suppressed in the presence of disorder. On the contrary, diagonal terms which contain orbit 
correlations on distances £, exhibit non-exponential behavior (Eq. (|5.7|) ) as a function of the 
inverse MFP l/l for integrable geometries and are not affected (within our approximations) 
by disorder in the chaotic case. 



C. Relation to experiment and other theories 

Measurements of the orbital magnetism of small microstructures are still rare today. The 
only experiment on ensembles of ballistic billiards that we are aware of, was performed by 
Levy et al.@ and investigated the magnetic susceptibility of an array of about 10 5 ballistic 
square-like cavities. The size of the squares is on average a = 4.5/^m, with a large dispersion 
(estimated between 10 and 30%) along the array. Each individual square is a mesoscopic 
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ballistic system since the phase-coherence length is estimated to be L$ = 15-40 \xm and the 
elastic mean-free-path I = 4.5-10 jum. The potential correlation length can be estimatedH 
to be of the order of £/a ~ 0.1. Taking the most unfavourable case of I ~ a ~ 4.5/im we 
obtain, with respect to the clean case, a disorder reduction for the averaged susceptibility of 
(x) fx — 0.37, showing that the features of the clean integrable systems (strong paramagnetic 
susceptibility at H = 0) persist upon inclusion of disorder. Since x — 100 Xl our 
calculations for the paramagnetic response of the ballistic squares agree quantitatively with 
the experimental findings (given the experimental uncertainties). 

Persistent currents in individual quasi-ballistic rings have recently been measured^. A 
similar setup would be needed for measuring the magnetic response of singly connected 
geometries, where our typical susceptibility ( |5.1|) should be measured for the integrable 
case. Since modern lithographic techniques allow to design chaotic as well as integrable 
cavities!!! and since we have demonstrated that disorder does not mask this difference, an 
order-of-magnitude effect is expected in the susceptibility according to the shape (chaotic 
vs. integrable) of the cavity. 

In a related theoretical work Gefen et al.il followed a complementary approach to ours 
and calculated the disorder-averaged susceptibility for an ensemble of ballistic squares based 
on long trajectories [strongly] affected by scattering from (5 like impurities. They found that 
the average susceptibility does not depend on the elastic MFP. These results are not borne 
out by either our analytical or our semiclassical calculations at temperatures relevant for the 
experiment, where the exponential damping from Eq. ( |2.18|) makes very long trajectories 
irrelevant. 

VI. SUMMARY 

In this work we have studied the interplay between integrability and disorder in the 
ballistic regime. The integrable property of the confining potential of a microstructure 
implies a peculiar behavior of its thermodynamical response functions, like the magnetic 



28 



susceptibility. The disorder effects provided by remote impurity scattering tend to weaken 
the importance of the boundary scattering (and therefore the relevance of the underlying 
classical mechanics). Using a semiclassical approach we quantify this damping and show it 
to be much weaker than previously estimated (power-law suppression instead of exponential 
damping for the typical and average susceptibility). The disorder damping is decisively 
affected by finite-size effects since it depends not only on bulk-like characteristics of the 
disorder (like the elastic mean- free-path), but also on the ratio between the size of the 
structure and the correlation length of the potential. 

Our finding for the weak disorder damping is particularly important due to the large 
phase coherence effects found for clean integrable structures and to the fact that the dif- 
ference in the magnetic response between integrable and chaotic geometries has not been 
experimentally demonstrated yet. 

Our calculational tools have been semiclassical expansions, which naturally convey at the 
quantum level the information about the underlying classical mechanics and its sensitivity 
with respect to disorder. For the weak disorder that we have considered in this work, the 
lowest order approximation consists of the perturbative modification of the classical actions 
by the impurity potential. Averages over impurity configurations following our semiclassical 
calculations, allow us to obtain various ensemble susceptibilities. Our analytical calculations 
have been checked against numerical quantum simulations performing exact diagonalizations 
of the corresponding Hamiltonian. 

The need to consider different averages is inherent to ballistic nanostructures, which 
are sufficiently small to be non-self-averaging. These various types of impurity-averaged 
susceptibilities for integrable systems are summarized in Table I for the three regimes defined 
by the correlation length of the impurity potential. We have first studied the fixed-size 
averaged susceptibility, directly obtainable from the disorder average of one-particle Green 
functions. It corresponds to the case where different impurity configurations of a given 
sample with a fixed Fermi energy are considered. For the short range regime, where the 
disorder correlation length £ < A F , we have an exponential suppression of the clean results 
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governed by the short-range elastic mean-free-path Is and the length of the most relevant 
trajectories. This result also holds in the finite-range (A F < ( « a), but with an elastic 
mean-free-path that we have evaluated semi-classically. In the long-range regime (£ > a) the 
fixed-size averaged susceptibility depends exponentially on the product (L / T) ■ (L / '£) (where 
L denotes the typical orbit length) and a correction taking into account the geometry of the 
periodic trajectories. 

For comparison with actual experiments we have to take into account that different 
impurity realizations are obtained together with a change in the Fermi energy and the size 
of the structures. We are then lead to consider impurity and size averaged susceptibilities, 
which are expressed in terms of two-particle Green functions. The typical susceptibility 
is appropriate when considering the magnetic response of an individual sample which is 
thermally cycled in order to obtain different realizations of the potential. The average 
susceptibility is obtained from the measurement of an array of microscopically different 
samples. For the short-range case the only difference between one- and two-point Green 
function quantities is the factor 1/2 of the exponential damping of the former. In the 
finite-range regime there appear important differences when considering two-point Green 
function quantities with respect to the one-particle case. Closed trajectories that remain 
closer than the correlation length of the potential result in a weak damping with a power- 
law dependence on l/L and £/a. This is the experimentally relevant situation, and the use 
of standard parameters lead us to conclude that disorder damping in currently realizable 
microstructures is sufficiently weak in order not to mask the large effects due to integrability. 
In the long-range case the damping due to disorder is extremely small. 

We have further considered the interplay between disorder and magnetic field in in- 
tegrable geometries. It is interesting to notice that both have a similar effect since they 
produce de-phasing between nearby trajectories. Since the two sources of de-phasing do not 
superpose, we find that disorder is less effective at finite fields, and reciprocally, disordered 
samples are less sensitive to magnetic field. 

In chaotic geometries periodic trajectories are usually isolated, resulting in smaller os- 
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cillations of the density of states and a much smaller magnetic response than integrable 
systems. Introduction of disorder in chaotic geometries is therefore less dramatic than in 
integrable systems, since it merely changes the action of the relevant periodic trajectories 
instead of producing de-phasing within a family. The transition from the ballistic regime 
(where classical trajectories are essentially unaffected by disorder) to the diffusive regime 
will be considered in a subsequent publication. 

In this work we have started from a system that is physically realizable using modern 
technology and we have developed a theoretical model with some key ingredients involving 
integrability and disorder. These are deep theoretical issues that need to be complemented by 
the consideration of other effects, like interactions, in order to obtain a complete description 
of the thermodynamics of mesoscopic systems. 
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APPENDIX A: RELATION BETWEEN SEMICLASSICAL AND QUANTUM 
MECHANICAL RESULTS FOR BULK MEAN FREE PATHS 



It is instructive to compare the semiclassical results of Eqs. ( |3.12|) -( |3.14|) for the ballistic 



regime with their counterparts obtained from quantum mechanical scattering theory. 

In a perturbative diagrammatic approach (treating the related Dyson-equation for scat- 
tering within a self-consistent Born approximation) the damping of the disorder-averaged 
one-particle Green function in a random potential is of the same exponential form as in 
Eq. (imp ! This is usually obtained by replacing the imaginary part of the self-energy in 
the Green function after impurity average by the density of states of the unperturbed system. 
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The resulting quantum mechanical inverse elastic MFP l qm , which appears in Eq. fl3.11|) , is 
related to the total cross section a by means of 

— = m(T, (Al) 

where n, is the impurity density and 

a = J d6 tr(6) (A2) 

with cr(G) being the partial cross section for scattering with an angle 0. 

For a Gaussian disorder potential of the form of Eq. (|3.3|) a calculation of the cross 
section can be performed analytically and the corresponding inverse MFP gives 

J- = - J„[2(A*) 2 ] e- 2 ^) 2 . (A3) 

'qm '<5 

Here, Iq is a modified Bessel function and 

1 2n nM 2 , niU 2 m ,. , 

l s h v F v F h A 

is the inverse MFP for the white noise case of 5— like scatterers of mean strength u. The v F 
is the Fermi velocity and <i(/i) = m/(27ih 2 ) the density of states at the Fermi energy of a 

2deg!. 

In order to compare l qm with our semiclassical result we expand l qm (k^) for large fc£ 
which gives 

i 

for k£ — > oo . (A5) 



l qm (kO ~ V47T (kO I 



The leading order term is exactly the semiclassical MFP Eq. ( |3.12|) for the Gaussian disorder 



model ( |3.3| ). The agreement between the semiclassical and diagrammatic approaches for the 
bulk can be related to the fact that our semiclassical treatment of disorder corresponds to 
the use of the Eikonal approximation (for each single scattering event) which is known to 
give the same results as Born approximation for large k£. 
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In the limit of £ < where our semiclassical description is no longer applicable, the 
mean free path l qm approaches lg, which means that Eq. ( |3.11| ) can further be used, but with 
the semiclassical I replaced by Is- 

The quantum mechanical transport mean free path It is calculated by including a factor 
(1 — cos 9) in the integral ( |A2[ ) for the scattering amplitude. It reads for Gaussian disorder 

i = \ {h\mf\ - h\m?\) e' m)2 (A6) 

for k£ — > oo . (A7) 



This relation shows that /t can be considerably larger than l qm for A F < £. This shows that 
in the case of a confined system and smooth disorder, the system may behave ballistically 
although the elastic MFP I might be considerably smaller than the system size. 
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FIGURES 

FIG. 1. Two representative periodic orbits characterized by xq and x' belonging to the family 
M = (1,1) (denoting, respectively, one bounce with each wall) of a square billiard of length a. 

FIG. 2. Magnetic susceptibility (x) (normalized with respect to the Landau susceptibility xl 
of a square billiard as a function of fcpa for the clean case (dotted) and for the ensemble average of 
billiards of fixed size with increasing Gaussian disorder (£/a = 0.1) according to an elastic mean 
free path l/a = 4,2,1,0.5 (solid lines in the order of decreasing amplitude). The susceptibility is 
calculated for zero magnetic field and at a temperature equal to 6 level spacings. 

FIG. 3. Logarithm of the ratio (x)/Xd as a function of the inverse elastic MFP a/l. The 
symbols indicate the numerical quantum results (from the top for £/a = 4, 2, 1, 0, 0.5 and 0.2. The 
dotted lines show the semiclassical analytical results for £/a = 4,2,1 (from above) according to 
Eq. ( p,16| ). The full line is the semiclassical result for £ = (Eq. ( |3.11 )). The quantum results 



for £ = 0.5 (squares) and 0.2 (diamonds) are beyond the regime of validity of the analytical limits 
f/a > 1 and £/a <C 1. 

FIG. 4. Averaged magnetic susceptibility (at H ~ of an ensemble of square billiards with 
variations in the size and impurity potential (£ = 0) for different disorder strength, i.e. elastic 
mean free path 1$. The full curves show the numerical quantum results and the dotted lines 
the semiclassical predictions from Eq. ( |5.6D taking into account the variations of lg with k F (see 
Eq. ( |A4[ )). The two sets of curves correspond to an elastic MFP lg/a = oo, 8, 4, 2, 1 (at k P a = 65), 
(from the top). 

FIG. 5. Logarithm of the ratio between disorder averaged and clean results for (a) typical 
X^ (b) ensemble averaged (x) susceptibilities as a function of increasing inverse elastic MFP a/l 
for different values of £/a. The symbols denote the numerical quantum results, the solid lines (for 
£ > 0) the semiclassical integrals (|5.l|) (a) and (|5.4|) (b) and the dashed lines asymptotic expansions 
(|5.9|) of the integrals for large a/l. 
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FIG. 6. Typical susceptibility as predicted by Eq. ( |5.1|) as a function of the dimensionless flux 
ip = Ha 2 /&Q. Dash line: clean case; solid line: / = a and £ = 0.1. 
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TABLES 





short-range 


finite-range 


long-range 


(x)/Xd 


exp (-L u /2ls) 


exp(-L n /2/) 


exp{-d 1 (L 2 //e)[l-/t/(2e 2 )]} 


(x^/x^) 2 


exp (-L u /l s ) 


ctit/aW/Lu) 1 / 2 


1 - d 2 a/l(a/0 9 


(x)/x 


exp {-Lu/ls) 


c a (Z/a)(l/L u ) 1/2 


1 - d 2 a/l(a/Z) 9 



TABLE I. Summary of the different average susceptibilities (at H = 0) considered in the 
short-range (£ < Xf < a), finite-range (Xf < £ < a) and long-range (Xf < a < £) regimes. The 
fixed-size impurity averaged susceptibility (x) is given by the one-particle Green function, while 
the typical an d average (x) susceptibilities are given by two-particle Green functions and 
involve impurity and energy averages. The different average susceptibilities are normalized with 
respect to the corresponding clean counterparts. L\\ is the length of the shortest flux-enclosing 
periodic trajectories in the square. In the short-range regime the damping is governed by the elastic 
mean- free-path lg given by the quantum mechanical expression ( |A4j ). The damping in the finite 
and long-range regimes is governed by the elastic MFP I (whose semiclassical expression is given 
in Eq. (|3.140 ), the correlation length £ of the impurity potential and the size a of the structure. 
It is the moment of inertia of the (11) trajectories (Eq. ( 3.1 6| ) ) . The finite-range expressions for 
X® and (x) showing a power-law damping hold in the deep ballistic limit I < L\\. The numerical 
factors are c t = (20/7)\/27r, c a = 2y/2^,d 1 = 1/40F, and d 2 = 6.5 • 1(T 5 . 
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Figure 1 




